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ON THE PERIODIC SOLUTIONS OF DISCONTINUOUS 
PIECEWISE DIFFERENTIAL SYSTEMS 

JAUME LLIBRE^ AND DOUGLAS D. NOVAES^’^ 


Abstract. Motivated by problems coming from different areas of the applied 
science we study the periodic solutions of the following differential system 

x'{t) = Fo{t,x) + eFi{t,x) + e‘^R{t,x,e), 

when Eq, Fi, and R are discontinuous piecewise functions, and e is a small 
parameter. It is assumed that the manifold Z of all periodic solntions of the 
unperturbed system x' = FQ(t,x) has dimension n or smaller then n. The 
averaging theory is one of the best tools to attack this problem. This theory 
is completely developed when Eq, Fi and R are continuous functions, and also 
when Fq = 0 for a class of discontinuous differential systems. Nevertheless does 
not exist the averaging theory for studying the periodic solutions of discontin¬ 
uous differential system when Fq ^ 0. In this paper we develop this theory for 
a big class of discontinuous differential systems. 

1. Introduction and statement of the main results 

1.1. Introduction. The study of the existence of invariant sets, particularly, 
periodic solutions is very important for understanding the dynamics of a differen¬ 
tial system. One of the most important tools to detect such sets is the averaging 
theory. A classical introduction to this tool can be found in [29l [28] . 

On the other hand the study of the discontinuous differential systems has it 
importance and motivation lying in some fields of the applied sciences. Many 
problems of physics, engineering, economics, and biology are modeled using dif¬ 
ferential equation with discontinuous right-hand side. For instance we may cite 
problems in control systems [2], impact and friction mechanics [1], nonlinear os¬ 
cillations mm , economics [El El, and biology mm- Recent reviews appeared 
in [301123]. 

Despite to the importance of the discontinuous differential systems mentioned 
above, there still exist only a few analytical techniques to study the invariant 
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sets of discontinuous differential systems. In [I9] the averaging theory has been 
extended for the following class of discontinuous differential systems 


( 1 ) 


x'{t) 


eFi{t, x) + X, e) if h{t,x)>0, 

eF 2 {t,x) + e‘^R 2 {t,x,e) if h{t,x)<0. 


where Fi, F 2 , i?i, R2 and h are continuous functions, locally Lipschitz in the vari¬ 
able X, T-periodic in the variable t, and h is a function having 0 as a regular 
value. The results stated in ra have been extensively used, see for instance the 
works [inillTlEHlEIlES]. 

In this paper we focus on the development and improvement of the averaging 
theory for studying periodic solutions of a much bigger class of discontinuous 
differential systems than in ([1]). Regarding to the averaging theory for hnding 
periodic solutions there are essentially three main theorems. In what follows we 
describe these theorems. 

The hrst one is concerning about the study of the periodic solutions of the 
periodic differential systems of the form 

x' = eFi{t,x) + e‘^F2{t,x) H-h e'^Fm{t,x) F e^^^R{t,x,e), 

with X G M'^. For continuous differential systems, even for the non-differentiable 
ones, this theory is already completely developed (see for instance [221 ESI El ESI 
[20]h and for discontinuous differential systems this theory is develop up to order 

2 in £ (see mm)- 

The other two theorems go back to the works of Malkin [23] and Roseau [27] . 
They studied the periodic solutions of the periodic differential systems of the 
form 


x' = Fait, x) + eFi(t, x) + e‘^F2{t, x) H-h e'^Fmit, x) + e^~^^R(t, x, e), 

with X G distinguishing when the manifold Z of all periodic solutions of 
the unperturbed system x' = FQ{t,x) has dimension d or smaller then d. These 
theories are well developed for continuous differential systems (see for instance 
[3 m [I2l ED])- Nevertheless there is no theory for studying such problems in 
discontinuous differential systems. Thus our main objective in this paper is to 
develop these last theorems for a big class of discontinuous differential systems. 

In subsections 11.21 and 11.31 we describe the class of discontinuous differential 
systems that we shall consider in this paper together with our main results, and 
we also do an application. In section [2] we prove our main results, and in section 
[ 3 ] we describe carefully the application of our results. 


1.2. Preliminaries. We take the ODE’s 

(2) x'{t) = F"'{t,x), (t,x)G§^xZl for n = 1, 2,..., A^, 
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where C is an open snbset and = M/T for some positive real nnmber T. 
Here x D — > for n = 1, 2,..., iV are continuous functions and the prime 

denotes derivative with respect to the time t. For n = 1, 2,..., iV let be open 
connected and disjoint subsets of x D. The boundary of S'„ for n = 1, 2,..., iV 
is assumed to be piecewise C™ embedded hypersurface with m > 1 and the union 
of all these boundaries is denoted by S. Moreover we assume that S and all 
Sn together cover x D. We call an N-Discontinuous Piecewise Differential 
System, or simply a DPDS, when the context is clear, the following differential 
system 


(3) 


F^{t,x) {t,x)^Si, 

F‘^{t,x) {t,x)eS 2 , 


x) {t, x) G S'at. 


Here Sn denotes the closure of in x D. 

Instead of working with system ([3]) we can work with the following associated 
system. 


N 

(4) x\t) = F{t, x) = ^ x^{t, x)F^{t, x), (t, x) X D, 

n=l 


where for a given subset H of x D the characteristic function XAif, x) is defined 
as 


XAif.x) 


1 if {t, x) G A, 
0 if (f, x) ^ A. 


Systems ([3]) and (jl]) does not coincides in S. Indeed system ([3]) is multivalued in 
S whereas system (|1]) is single valued in S. Using Filippov’s convention for the 
solutions of the systems ([3]) or (|1]) (see [n] ) passing through a point (t, x) G E 
we have that these solutions do not depend on the value F(t, x). So the solutions 
of systems ([3]) and dH) are the same. 

When F^ for n = 1, 2,..., iV are functions we dehne the “derivative” of the 
discontinuous piecewise differentiable function F{t,x) with respect to x as 


N 

(5) D^F{t, x) = ^ x)D^F^{t, x). 

n=l 


We note that when the function F{t,x) is differentiable with respect to the vari¬ 
able X then the above definition coincides with the usual derivative. 


We say that a point p G S is a generic point of discontinuity if there exists a 
neighborhood Gp of p such that iSp = Gp n S is a C™' embedded hypersurface in 
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E)^xD with m > 1, such that the hypersurface Sp splits Gp\Sp in two disconnected 
regions, namely Gp and G~, and the vector helds F+ = F\q+ and F~ = F|g- 
are continuous. We dehne l{p) as the segment connecting the vectors Fp{p) and 
F~{p) when these have the same origin (see Figures [Hand [2]). 

Let iS C S be an embedded hypersurface in x D and TpS denotes the tangent 
space of S at the point p. In what follows we dehne the crossing region 
(see Figured]), and the sliding region S'^(iS) (see Figure |2|) of the hypersurfaced. 

S"(5) = {peS: l{p) n TpS = 0} and E"(5) = {p e S : l{p) n TpS ^ 0} . 

These dehnitions only make sense when the linear space TpS is based at the origin 
of the vectors Fp{p) and F~{p). 

The hypersurface iS C E can be decomposed as the union of the closure of its 
crossing region with its sliding region. 

When the hypersurface iS C E is given by iS = h“^(0) for some function 
h : X D —!■ M having 0 as a regular value, then the above dehnitions becomes 

E'(S) = {p e 5 : (V/!(p), (1, F*(p))){Vh(p), (1, F-{p))) >0} and 
E‘(5) = {p e 5 : {Vft(p), (1. F+(p))>(V;i(p). (1, F-(p))) < 0} . 




Figure 1. Crossing region of <S: E'^5. 

Globally we dehne the crossing region YT as the generic points of discontinuity 
p such that p G E‘^(iSp). The sliding region E^ is dehned analogously. Later on 
this paper for a point p G E'^ we shall denote TpE = TpSp. 

Let ipF^{t,q) be the solution of system ([2|) passing through the point q E Sn 
at time f = 0, i.e. <pipn(0,g) = q. The local solution ippit^q) of system (|1|) 
passing through a point p G E'^ at time f = 0 is given by the Filippov convention, 
i.e. for p G E'^ such that l{p) C G^ and taking the origin of time at p, the 
trajectory through p is dehned as ppit^p) = Tp-ii^yP) t E Ip H {t < 0}, and 
Lfpit^p) = ipp+{t,p) for t E Ip 0 {t > 0}. Here Ip is an open interval having the 
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Figure 2. Sliding region of S: 

0 in its interior. For the case l(p) C G~ the dehnition is the same reversing the 
time. 

Assuming that the functions F"'{t,x) are Lipschitz in the variable x for n = 
1, 2,..., iV, the results on Filippov systems (see m) guarantee the uniqueness 
of the solutions reaching the set of discontinuity only at points of 

1.3. Statements of the main results. Let D be an open subset of and 
for n = 1, 2,..., iV let x D —)■ be a C™' function with m > 1, and 

F” : X H —)■ and x D x [0,1] —?■ be continuous functions which 

are Lipschitz in the second variable. All these functions can be seen as T-periodic 
functions in the variable t when f G R. Later on in this paper we shall assume 
more conditions under these functions. 

Now taking 

N 

Fi{t, x) = ^ x^{t, x)Fp{t, x), for f = 0,1, and 

71=1 

N 

R{t,x,e) = 

n=l 

we consider the following DPDS, 

(6) x'{t) = Fo(f, x) + sFiit, x) + e^R{t, x, e). 

The parameter e is assumed to be small. We recall that E denotes the union of 
the boundaries of Sn for n = 1, 2,..., iV. 

In order to present our main results we have to introduce more dehnitions and 
notation. 

For z ^ D and e > 0 sufficiently small we denote by x{-,z,e) : [0,^(2^^)) —)■ R'^ 
the solution of system ([6]) such that x{Q,z,e) = 2 ;. Given a subset F of F we 
dehne = {(f, x(f, z, e)) : z & B,t E [0,f(z,e))}- 
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We denote by Sq the set of points x E D such that the function F{0,x) is 
discontinuous, clearly {0} x Sq C E. 

One of the main hypothesis of this paper is that the unperturbed system 

(7) x'{t) = Fo{t,x), 

has a manifold Z embedded in D\d'EQ such that the solutions starting in Z are 
all T-periodic functions and reach the set of discontinuity S only at its crossing 
region Here (9Eo denotes the boundary of Eq with respect to topology of D. 
Precisely, 

(H) let Z = {za = {a, Po{a)) : a G H}, where V is an open and bounded 
subset of and (3o : V ^ is a C"* function with m > 1. We shall 

assume that Z G D, Z H (9Eo = 0, O E C E'^ and for each z^E Z the 
unique solution Xa(t) = x(t,Za,0) is T-periodic. 

Remark 1. Suppose that the solution Xa{t) reaches the set E^ times. The 
assumption Z n c^Eq = ^ in hypothesis {H) implies that for each z^ E Z there 
exists a small neighborhood Ua G D of Za such that for e > 0 sufficiently small 
every solution of the perturbed system ([6]) starting in Ua reach the crossing region 
of the set of discontinuity E^ also Aq times. This fact will be well justified in the 
proofs of Lemmas\B and\^ in sectionl^ 

For z E D we take the following discontinuous piecewise linear differential 
system 

(8) y'= D^Fo{t,x{t,z,0))y, 

which can be seen as the linearization of the unperturbed system ([7]) along the 
solution x(t, z, 0). We note that for each z E D the function 1 1 —>■ D^Foit, x(t, z, 0)) 
is piecewise C™' with m > 1, so we can consider a fundamental matrix Y [t, z) of the 
differential system ([H]). Clearly t i—)■ Y{t,z) is continuous piecewise C"* function. 
We dehne 

(9) yi{t,z) =Y{t,z) [ Y{s,z)~^Fi{s,x{s,z,0))ds. 

Jo 

Now for Za E Z we denote Yaifi ) = Y { t , Za ). Let vr : x R'^“^ —)■ R^ and 

7r-*“ : R^ X R"^”^ —)■ R"^”^ be the projections onto the hrst k coordinates and onto 
the last d — k coordinates, respectively. Thus we dehne the averaged function 
/i : F ^ R^ as 

(10) fi{a) = 7ryi{T,Za). 

In what follows dis(x. A) denotes the Hausdorff distance function between a 
point X E D and a set H C T, and as usual the function dB^fi^W^tf) denotes 
the Brouwer degree (see for instance [5] for details on the Brouwer degree). Our 
main result on the periodic solutions of DPDS ([6]) is the following. 
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Theorem A. In addition to the hypothesis {H) we assume that 

{HI) forn = 1,2,... ,iV, the functions Fq and (Iq are of class C^; the continuous 
functions D^Fq, Ff and R are locally Lipschitz with respect to x; and the 
boundary of Sn are piecewise embedded hypersurface inM. x D; 

{H2) for any a E V there exists a fundamental matrix solution Y{t, z) of (jH]) 
such that the matrix Ya{T)Ya{0)~^ — Id has in the upper right corner the 
null kx {d—k) matrix, and in the lower right corner has the {n—k) x {n—k) 
matrix A„ with det(AQ,) ^ 0 ; 

{H3) for an open subset U of D such that Z <Z U we have that (0, i/i(s, 2 ;)) G 
T{s,x{s,zfi))'^ whenever (s, a:(s, z, 0)) G for (s, z) G x U; 

{HY) there exists W open subset of V such that fi{a) 7 ^ 0 for a G dW and 
dB{fi,W,0)^0. 


Then for e > 0 sufficiently small, there exists a T-periodic solution (p{t,e) of 
system ([H]) such that dis(<y5(0, e), Z) ^0 as e —)■ 0. 


Theorem is proved in Section [2l 


Remark 2. When fi is a W function the assumption 

(h4) there exists a eV such that fi{a) = 0 and det(/((a)) 7 ^ 0, 
is a sufficient condition to guarantees the validity of the hypothesis {HA). 


Theorem B. We suppose that the hypotheses {H), {H2) and {H3) of Theorem 
HI hold. If we assume that 

{hi) for n = 1,2,..., N, Fq, D^Fq, Ff, , and fo are functions and the 

boundary of Sn are piecewise embedded hypersurface inM x D, 

then fi{a) is a W function for every a eV. Moreover, if we assume in addition 
that hypothesis (M) holds, then for e > 0 suffciently small, there exists a T- 
penodic solution (p{t,e) of system ([H]) such that (p{ll,e) Za as e —)■ 0 . 


In what follows we provide an application of Theorems and [Bl We study 
the existence of limit cycles which bifurcate from the periodic solutions of the 
linear differential system {u,v,w) = {—v,u,w) perturbed inside the class of all 
discontinuous piecewise linear differential systems with two zones separated by 
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( 11 ) 


V 

w 


= 


\W J 


/ 


the plane S = {n = 0} C i.e. 

^ —V + e{af + bfu + cfv + dfw) ^ 

u + e{a2 + b^u + C 2 V + d^w) if n > 0, 

y w + e{a^ + b^u + c^v + d^w) j 

—V + e{ai + b^u + c{v + d^w) ^ 

u + e(a^ + b^u + c^v + d^w) if n < 0. 

y w + e{a^ + b^u + c^v + d^w) y 

Our result on the existence of a limit cycle of system ffTTj) is the following. 

Proposition 1. If (af — a^)(6h + ^ 1 " + cf + cf) > 0, then for |e:| > 0 sufficiently 
small there exists aperiodic solution {u(t,e),v(t,e),w(t,e)) of system flTT]) such 
that ta(0,£) 0 when e —?■ 0. Moreover, we can find {u*,v*) G such that 

7r(6i +b{ +C 2 +cf) 

and (m( 0, e), n(0, e)) —)■ {u*, v*) when e —)■ 0. 

Proposition [U is proved in Section |3l 


2. Proof of Theorem 


Before proving our main result we state some preliminary lemmas. 

Given a function f : [0,1] we say that f{e) = O(e^) for some positive 

integer £ if there exists constants £1 > 0 and k > 0 such that ||^(e)|| < k\e^\ for 
0 < £ < £ 1 , and that .^(e) = o{e^) for some positive integer £ if 


lim 

£—^0 


lle(g)ll 


0 . 


Here 11 ■ 11 denotes the usual Euclidean norm of R'^. 
called the Landau’s symbols (see for instance [28]). 


The symbols O and o are 


Lemma 2. Under the hypotheses {H), {HI), and {H3) of Theorem\M there exist 
an open and bounded subset C of U\dTiQ, a compact subset Z (Z C with Z C Z°, 
and a small parameter > 0 such that t(^z,e) > T for every z E C and e G [0, £ 9 ] • 
Moreover x{t, z,e) = x{t,z,0) +eyi{t,z) + o(e) for every z E Z, e E [0,eo]; 
t E [0,T]. Here Z° denotes the interior of the set Z with respect to the topology 
of D, and the function yi is given in ([9]). 


Proof. We note that Z is a compact subset of D and SEq is a closed subset of D, 
such that, from the hypothesis {H), Z fl SEq = 0. So there exists an open subset 
A of D such that Z (Z A and A fl OEq = 0. 
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Also from hypothesis {H) we have that for a E V the continuous function Xa{t) 
reaches the set S only at points of Since this function is T-periodic we can 
hnd a hnite sequence (t^) for i = 0,1,..., with = 0 and = T such that 

xl{t) 


if 0 = 

tl<t< tl, 


if 


Xait) = < 


X^t) if 


x^-it) if 


< f < f 


ty ^ <t<ty 


= T, 


where each curve t i—>■ = x^(t, Za, 0) reaches the set only at t = and 

t = for i = 2, 3,..., Aq, — 1, the curve reaches the set only at t = 0 
and t = if (0, Za) E S, and only a.t t = if (0, z^) ^ S, the curve 
reaches the set only at t = and t = T if (T, x(T, Zq,, 0)) E S, and only 

at t = if (T, x(T, Zq,, 0)) ^ S. From the dehnition of the crossing region E^ 
these intersections are transversely. 

Since for i = 1, 2,..., are solutions of Lipschitz differential equations we 
can use the results of continuous dependence of the solutions on initial conditions 
and parameters to ensure the existence of a small parameter Sa and a small 
neighborhood Cq C Ant/ of Zq such that Ca HE C E^ for every e E [0, £q]. From 
the compactness of Z we can choose ei as a minimum element of E V. Now 
taking C = it follows that n E C E'^ for every e e [0,£i]. Moreover, 

we can take > 0 and C smaller in order that the function t i—)■ x{t, z, e) is 
dehned for all {t,z,e) E x C x [0,£i]. This is again a simple consequence of 
the Theorem of continuous dependence on initial conditions and parameters. 

Thus for 2 ; G C and £ E [0,ei] the function t i-A x{t,z,e) is continuous and 
piecewise Ch So we can hnd a hnite sequence {f{z,e)) for i = 0,1,... with 
t^{z,e) = 0 and t'^^{z,e) = T such that 


x^(t, z, e) 
x‘^(t, z, e) 


if 0 = 
if 


(12) x(t, z^e) = < 


x^{t,z,e) if 


£) Zt < t^{z, e), 
s) <t < t^{z, e), 

f~^{z,e) <t< f{z,e), 


^{z,e) <t< = T, 


x‘^^{t,z,e) if 
for which we have the following recurrence 

(13) x^{0,z,e)=z and x\f~^{z,e), z,e) = x’‘~^{f~^{z,e), z,e), 
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for i = 2, 3,, Kz- The crossing region is an open subset of S, so for each 
z ^ C we can hnd 0 < < Si sufficiently small such that the number of 

intersections between the curve t i-A- x{t, z, e) with the set S'^ for 0 < t < T and 
for e G [0, Ez] does not depend of e. Since C is compact we can find 62 a minimum 
element of the e^’s for z & C such that the above statement holds for every z & C 
and e G [0, £ 2 ]- 

Here again for every z E C and £ G [ 0 , 62 ] each curve t i-A- x\t,z,e) reaches 
the set E^ only at t = f~^{z, e) and t = e) for i = 2, 3,..., — 1, the curve 

z,e) reaches the set E^^ only at t = 0 and t = t^{z,e) if (0, z) G E, and 
only at t = t^{z,e) if ( 0 , 2 ;) ^ E, the curve x^^{t,z,e) reaches the set E^ only 
at t = t'^^~^{z,e) and t = T if (T, x(T, z, 0)) G E, and only at t = t'^^~^{z,e) if 
(T,a:(r,z,0)) i E 

The functions t 1 —)■ x^{t, z, e) for i = 1,2,..., are and satisfy the DPDS 
dHD, so there exists a subsequence (uj) for z = 1,..., with n* G {1,2,..., N} 
such that 


(14) ^x\t,z,e) = FQ'{t,x\t,z,e)) + eF^^{t,x’‘{t, z,e)) + e'^R"‘^{t,x\t, z,e),e). 


Therefore the function x\t,z,e) is the solution of the Cauchy Problem defined 
by the differential system (114)) together with the corresponding initial condition 
given in flTS]) . Moreover x*(t, Za, 0) = and f{za, 0) = for z = 1, 2,..., 

From the continuity of the function x{t, z, e) we can choose a compact subset 
K oi D such that x{t,z,e) G K for all {t,z,e) G x C x [0,£2]- From the 
continuity of the functions F” and for z = 0,1 and n = 1,2,..., N we have 
that these functions are bounded on the compact set x K x [0, £ 2 ], so let M be 
an upper bound for all these functions. Let L be being the maximum Lipschitz 
constant of the functions F", DFq, and F” for z = 0,1 and n = 1,2,... ,N on 
the compact set x F x [ 0 ,£ 2 ]- 

We compute 


R{s, x(s, z, e),e)ds 


< 


I|F(s, x{s, z,e),e)\\ds = TM, 


R{s,x{s, z,e),e)ds = 0{1) in the parameter £. 


which implies that 


0 
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For z E C and t G (0, T) we can find k G {1, 2,..., — 1} such that t G 

e),£)) and 


x{t,z,e)= x'^{t,z,e) 


= x'" ^{z,e),z,e) + 




Fo(s, x(s, z, £:))ds 






Fi{s, x(s, z, e))ds + 0 {e‘^). 


Since 


x^{f{z,e),z,e)= X* e),z,e)+ 


rt*{z,e) 


-\-£ 


p-^{z,e) 

Fi{t, x(t, z, £))(it + 0 {e‘^), 


Fo(t, x(t, 2 ;, e))dt 


for i = 1, 2,..., we obtain, proceeding by induction on i, that 

(15) x{t,z,e)=z + / Fo{s,x{s, z,e))ds + e / Fi{s,x{s, z,e))ds + 0{e^). 


From here the proof of the lemma follows by proving several claims. 


Claim 1. There exists a small parameter £0 > 0 such that for any z ^ Z and for 
i = 0,1,2,Kz the function f{z, e) is of class for every z in a neighborhood 
Uz C C of z and for e G [0,£o]; {df/de){z,0) = 0. Moreover for T~^{z,0) < 

t < t\z, 0) we have that yi{t, z) = {d x'' /de){t, z, 0) for i = 1,2 ,..., Kz. 

First of all we note that F{z,e) = 0 and F^{z,e) = T. So the hrst part of 
Claim [T] is clearly true for i = 0 and i = Kz- 

We have concluded above that for each z & Z the curve t 1 —)■ x{t, z, 0) reaches 
the discontinuity set only at points of Let z^ = x\f{z,0),z,0) and p* = 
{f{z, 0), z^) G E'^, then p* G E'^ for every i = 1,2,..., Kz il (0, x(T, 2 ;, 0)) G E, and 
for every z = 1, 2,..., — 1 if (0, x(T, 2,0)) ^ E. Particularly pi is a generic point 

of E, so there exists a neighborhood Gpi^ of p* such that iSpi = Gp* n E is a 
embedded hypersurface of x H with m > 1. It is well known that iSp* can be 
locally described as the inverse image of a regular value of a function. Thus 
there exists a small neighborhood Gpi of pi with Gpi C Gpi and a C™ function 
hi : Gpi^ —)■ M such that Gpi n iSp* = hL^(O) fl E. 

For (t, x) G Gpi^ system ([6]) can be written as the autonomous system 

X{T,x,e) if hi{T,x)>0, 

Y{T,x,e) if hj(r, x)<0. 
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where 


X{r,x,e) -- 
Y{T,x,e) 


(r, x)eF^^+^ (r, x) + (r, x, e) 

1 

Fq”* (r, x) + eFf' (r, x) + (r, x, e) 




From the definition of crossing region we also have Xhi{pl,0)Yhi{pl,0) > 0, 
then 


(16) 


0 jt Yhi(p\, 0) = ^ < (1 A”'(A)) 


Now defining Hi(tXx) = hi{t,x''{tXx)) have that Hi(t''{z,0), z,0) 
Since 


= 0 . 


dH, 


d 


0), z, 0) = F(t, C, e)) 

dh. 




^(f (z,0),F(f(z,0),z,0)) 
f)h ■ 

0). 0), z, 0))^(i'(j, 0), 0) 




dh, 
dt 

t(pa + t(p;)A"‘(pa 


= r/!.(p‘,o)#o, 


from the Implicit Fnnction Theorem we conclnde that there exist a small neigh¬ 
borhood Uz <Z C oi z and a small parameter > 0 snch that e) is the nniqne 
C™ fnnction with = 0 for every C, eUz and e G [0,ez]. So 


f) 

(17) = F(C, 0) -I- 0) -I- 0(e) 

for every i = 1, 2,..., — 1. Since Z is compact we can take eo as a minimnm 

element of e^’s for z E Z. 
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Now we shall use finite induction to conclude the proof of Claim [H We note 
that hi{V-{z, e), z, e)) = 0 ioi e ^ [0, ^o]) so 


0 = £),x\f{z, e), z, e)) 


£=0 


dh, 


dh 


= o7(P^) —0) + 0), 0) — {z, 0) 


dx^, ,■, 


df 


dt 


de 


(18) 


dt de 


3h 3 3h { 33 3 t^ 

:b^) —0) +-(pl) ( Fo-(p;) —(z, 0) + —(f (z, 0), z, 0) 


dt ^ de 


dz 


/ / 3 3 -33 3t^ 

= (vh{pl), (—(z, 0), F^^{pl) — {z, 0) + —{f{z, 0), z, 0) 

for i = 1,2,, Kz. 

Taking i = 1, from flTT|) we obtain that 
(19) ^ ^ DFQ^{t,x^{t,z,ll)) 

So for 0 < f < f^( 2 ,0) the differential system (II Dll becomes 


(20) 17 


^ (^,^,0)^ = DFQ{t,x{t,z,Q)) (^^{t,z,0) ] + Fi{t,x{t,z,0)). 


dt V de 


Since -^—(0, 2 ;, 0) = 0 the solution of the linear differential system fl20|) is 


( 21 ) 


dx^ 

de 


{t,z,0) = Y{t,z) / Y{s,z) ^Fi{x{s,z,0))ds = yi{t,z). 


for 0 < t < t^{z, 0). Now from hypothesis {H3) and from equality fflSj) . for i = 1, 
we have that 

( 22 ) (^X—{z,0),XF-^{pl) — {z,0)+y,{t\z,0),z))eT,.J: 

for every A G [0,1]. Thus 

(23) 

0 = /v/r(p)), (a^(z,0),AF„"‘(p()^(z, 0) + ai(t'(z, 0), z)) \ 

/ f)h 3 3 3h 33 \ 3h 

3 3 3h 

= XYhipl, 0) — {z, 0) + —{pl)y,{t\z, 0),z), 
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for every A G [0,1]. Computing the derivative with respect to A in (l23i) it follows 
that Yhi{pl,0)-^{z,0) = 0 . So from flTdl) we obtain that 


(24) 


de 


{z,0) = 0. 


Hence from fl2Tll and fl2411 the claim is proved for i = 1. 

Given a positive integer i > 1, we assume by induction hypothesis that Claim 
[T] is true for i = i — 1. Taking i = i from flTT)) we have that 

(25) ^ +^r(^,a;^(C^,0)). 

So for t^~^{z, 0) <t < t^{z, 0) the differential system fl25|) becomes 


(26) 


d (dx^ 


dx^ 


dt = DFo{t,x{t,z,0)) (t,2r,0)J + Fi{t, x{t, z,0)). 

From flT^ we have that x^{t^~^{z,e), z,e) = x^~^{t^~^{z,e), z,e) for every e G 
[0,eo]- Computing its derivative with respect to e at e = 0 we obtain that 


(27) 


0 ), 0)^{z, 0 ) + 0 ), 0 ) = 

(t^-^(.,0),.,0)^(.,0) + ^^ 


(t^-^(z,0),z,0). 


dt de de 

So from induction hypothesis it follows that 

(28) _(t^-i(^,o),z,0) = -^{t^-\z,Q),z,Q)=y,{t^-\z). 

We note that (l28|) is the initial condition for system (126|) . Thus for 0) < 

t < t^{z, 0) regarding to the linear differential equation fl26ll we get that 


(29) 


de 


{t,z,0) = Y{t,z)yi{t^ ^{z,0),z) + Y{t,z) I Y{s,z) ^Fi{x{s, z,0))ds, 


where Y{t, z) is the fundamental matrix of the linear differential system ()5|) 
such that Y{t^~^{z,0), z) is the identity matrix. Clearly Y{t,z) = Y{t, z)Y{t^~^ 
( 2 , 0),^)“^. So substituting ([9]) in fl2^ we get 

dx^ 

— (t,z,0)= Y(t,z) Y{s,z)~^Fi{x{s,z,0))ds 


( 30 ) 


+Y{t,z) / Y{s,z) ^Fi{x{s, z,0))ds 

J t^“l(2,0) 

= Y{t,z) [ Y{s,z)~^Fi{x{s,z,0))ds = yi{t,z), 

7t^-i(2,o) 
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for <t< 0). Now repeating the procedure of (1^ and (j23|) for 

dt^ 

i = i we conclude that 0) = 0. So we have proved Claim [1] 

oe 


Claim 2. The equality x{t,z,e) = x{t,z,0) + 0{e) holds for every z ^ Z and 
e G [0, £o]- 

If t G [t'^~^{z,e),T^{z,e)) then we compute 
rt 

/ Fo{s,x{s, z,£))ds = / Fq^{s,x{s, z,£))ds 

Jq \Jp-^{z,e) 

+ f Fq'^{s,x{s, z,e))ds 

J 

«-i / pdizfi) 

= ^ / FQ^{s,x{s,z,e))ds 

i=i Vt*-Hz,o) 

+ [_ FQ^{s,x{s,z,e))ds + Eo{e), 


where 


'^-1 / i-r-^{z,o) 


rPizfi) 


Eo{£) = y] I / F^^{s,x{s,z,e))ds- I Fq^{s,x{s, z,e))ds 

i=l dp{z,e) 






Fq'"{s, x(s, z, £))ds. 


It is easy to see that there exists a constant E such that 


K—1 


(31) 




i=0 


Indeed the function Fq" (t, x) are bounded in the set x iC, so 


r-T( 2 , 0 ) 


' P{z,e) 


Fq"{s, x{s, z, e))ds 


< / \\F^'{s,x{s,z,e))\\ds 

J P[z,e) 

< L\t\z,t)) -t\z,e)\, 


for i = 0,1, 2,, K. 
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From Claim [U we conclude that Eo{e) = o{e), particularly Eq^e) = 0{e). Thus 


1 


Fo(s, x(s, z, e))ds = 


i>t\z,0) 


(32) 


i=l 


Eq"{s, x{s, z, e))ds 




x{s, z, e))ds + 0{e). 


Using the fact that the functions Eq" for i = 1,2,..., Kz are locally Lipschitz 
in the second variable together with (l3^ we obtain 


'0 

''“l rt'(z,0) 

X-l(2,0) 

I 

K-1 

^ / 

i=l dti-\z,0) 


-\-L 


Eq{s, x{s, z, e)) — Eo{s, x{s, z, 0))ds 

||Fo”'(s,a:(s,z,£)) - Fo”'(s, a:(s, 0)) | |ds 

Eq^{s, x{s, z, e)) - Eq^{s, x{s, z, 0))| |ds + 0{e) 
x{s,z,e) — a;(s, 2 ;, 0)||(is 
x{s, z, e) — x{s, z, 0)1 |ds + 0{e) 




= L \\x{s, z,e) — x{s, z,ll)\\ds + 0{e). 

Jo 

From flTHll this implies that 


(33) 


\x{t,z,E) - x{t,z,0)\\ < / ||Fo(s,a;(s,2;,£)) - Fo(s,x(s,2r,0))||cis 

Jo ^ 

+ kl f \\Ei{s,x{s, z,E))\\ds + 0 {e^) 


< I^IMT + L / I |x(s, 2 ;, e) — a:(s, 2 ;, 0)1 Ids 


< |e|MTe^^. 

The last inequality is a consequence of Gronwall Lemma (see, for example, Lemma 
1.3.1 of [28]). 

From fl5^ we conclude that x{t,z,E) = x{t,z,0) + 0 {e). So we have proved 

Claim [21 
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Claim 3. The equality x(t, z, e) = x(t, z, 0)+eyi(t, z) + o{e) holds for every z E Z 
and e G [0, ^o]- 

In the proof of Lemma 1 of [20] it has been proved that 


Fff\t,x\t,z,e)) = Fff\t,x\t,z,t))) + D^FQ\t,x\t,z,t))) 

(34) ■{x''{t, z, e) - x\t, z, 0)) + 0{e^), 

Ff*{t,x\t,z,e)) = Ff\t,x\t,z,t))) + 0{e), 

for all T~^{z, s) < t < T{z, e) and for every i = 1, 2,..., In what follows we 
give a sketch of the proof. 

Let >C(/i) = G{t, fix^{t, z, e) + (1 — /i)a:*(t, z, 0)). So 

(35) 

Fo"'(t, z, s)) = F^^ (f, z, 0)) + £i(l) - A(0) 

= F^'{t,F{t,z,0))+ ['c[{X,)dX, 


= F^^{t,F{t,z,0)) + 


D,F^^{t,G{F{t,z,e)))dXi 


{x'^{t, z, e) — x*(t, z, 0)) 


D,F^^ {t, e,{F{t, z, e))) - F{t, z, 0)) 


dAi 


Jo 

■{x\t,z,e) -x'‘{t,z,0)) + FQ^{t,x\t,z,0)) 
+D^Fo"{t,x\t,z,0)) ■ (x\t,z,s) - x\t,z,0)). 


So observing that the fnnction D^FQ*{t,x) is locally Lipschitz in the variable x 
and using Claim [2] in (1351) we obtain the equality for Fq^ of (1341) . The equality 
for Ff^(t,x) of (134|) is obtained directly by using Claim 12] together with the fact 
that this function is Lipschitz in the variable x. 

From (134|) we obtain that 


F^^{t,x\t,z,e)) = 


Fq* {t, x^{t, 2 ;, 0)) + eD^Fq* (t, x^{t, z, 0)) 

—(i,z,0) + C>(P), 


(36) 
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for all t* < t < t'^{z, e) and for every i = 1, 2,..., For the moment we 

dx^ 

cannot use Claim [1] to ensure that -^{t,z,0) = yi(t,z) because it is only true 
when f~^{z, 0) <t < ^{z, 0). 

Given z G C we have that, for every f~^{z,e) < t < f{z,e), x^{t,z,e) = 
x{t, z, e) for z = 1, 2,, Ka- Moreover if f~^{z, s) < s < t^{z, s) and e G [0, ^o]; 
then Fp{s, x\s, z, e)) = Fj{s, x{t, z, e)) for j = 0,1 and for every z = 1, 2,..., F. 

If F~^{z,e) <t< F{z,e) from fl3il) we compute 
(37) 

/ Fi{s,x{s, z,e))ds = 


«-l rt'{z,e) 


E 


=1 >zn-i(2,£) 


F”*(s, a:*(s, z, £))(is I + / F^'^{s,x'^{s, z,e))ds = 

J t^~^{z,£) 


' K.— 1 


-zy2,£) 


/ F”'(s, a;*(s, 2 ;, 0))(is j + / F^'^{s, x'^{s, z,0))ds + 0{e) = 

_i=l j Ji^-^{z,e) 

rt^iz,0) \ pt 

/ F"*(s, x*(s, z, 0))ds I + / Fl^'^{s,x'^{s, z,0))ds + Ei{e) 


i=i Jti-^{z,o) 

+0(e) = 


tK-l(^Z,0) 


/ Fi(s,a:(s,.s,0))ds ) + 

Jn-i(2,o) j 

ft 

Fi{s,x{s, z, 0))ds + Ei{e) + 0{e), 




Fi{s, x(s, z, 0))ds + Ei(s) + 0(e) = 


where 


Bi(£)= E 


K-1 / rt‘-^(z,0) 


rtdz,0) 


F^"(s, x’'(s, z, 0))ds — / F”'(s, a;*(s, 2 ;, 0))(is 

\yJt^-^(z,£) Jt^(z,e) 

|■i^-^{z,0) 

+ / F"^(s, z, 0))(is. 

J t^~^(z,e) 

Now, as in the case Eq^e) of the proof of Claim [21 it is easy to see that there 
exists a constant E such that 

K—1 




i=0 


(38) 
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So from Claim [U we conclude that Ei{e) = o(e) and consenquently Ei{e) = 0{e). 
Going back to inequality fn7|) we obtain 


(39) 


Ei{s^x{s^ z^e))ds = / Fi(s,a:(s, 2;, 0))(is + (9(e). 


dx^ 

From Claim [H ——(t, .^,0) = yiit^z) for f~^(z,0) <t< f(z,0), so from 
oe 

we compute 
(40)^ 

/ Fa{s,x{s,z,e))ds ^ 






=1 


Eq"{s,x’'{s, z,e))ds \ + / Eq'^{s,x^{s, z,e))ds = 




rt\z,e) 


K.— 1 

i=i Ft) 

rt 


dx" 


[F()*'(s,x*(s,z, 0))+eF>^Fo"'(s,x*(s,2;,0)) —(s,2;,0)]cis ) + 

dx^ 


K-l / |.T(2,0) 

S \Jti-Hz,o) 

rt 


[Eq'‘{s, x'"{s, z, 0)) + sD^E^'^is, x^(s, 2 ;, 0))-^(^> 0)]^^ + 0{e^) = 


dx'' 


[Eq*{s,x\s,z, 0)) + eD^Eo^{s,x\s,z,0)) — {t,z,0)]ds ) + 


dx'^ 


[Eq'^{s, x'"(s, z, 0)) + eD^E^'^is, a:'"(s, z, 0)) —(t, z, 0)]ds + F;2(e) 




Jtf^-^(z,o) ■ 9e 

+0(e^) = 

^ / /-PC.o) \ 

[Eo{s,x{s,z, 0 )) + eD^EQ*{s,x{s,z,0))yi{s,z)]ds + 

\Jy-^{z,o) J 

f [Fo(s,x(s,z, 0)) + eD^EQ'^{s,x{s,z,0))yi{s,z)]ds + E 2 {e) + 0(e^) 

The last equality comes from observing that Eq^{s,x\s, z,0)) = Eq{s,x{s, z,0)) 
for every s G 0), C(.s, 0)) and i = 1,2, From definition (jS]) the 

inequality fHOl) becomes 


( 41 ) 


Eo{s,x{s,z,e))ds = / [Fo(s, x(s, 2;, 0)) + eF>,i,Fo(s, a;(s, 2;, 0)) 

Jo 

'2/1 (s, z)\ds + E 2 {£) + (9(e^). 
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Here 


K,— l / 1 


B2(e)= Y. 




= 1 

A(2,0) 


[Fj*'(s, x*(s, z, 0)) + eD^Fq^s, x\s, z, 0))-^(^> 


(9a;* 


+ 


h^(z,£) 
/•t^-l(2,0) 

I t^~^{z,£) 


[Fo"'(s, a;(s, ^r, 0)) + sD^F^'is, x\s, z, 0))-^(^> 0)] 


9x* 


[Fo**"(s, x*(s, z, 0)) + £D^Fq'‘{s, a;*(s, z, 0))-^(^) 


It is easy to see that there exists a constant E such that 


K—1 


(42) 


I|B2(£)II<-B5^|«‘(2:,0)-A2:,£)|, 


i=0 


So from Claim [T] it follows that E 2 {£) = o{£). Going back to inequality fHT]) we 
have 


Fo{s, x{s, z, £))ds = / Fo(s, x(s, z, 0))(is 


(43) 


+£ / D^Fo{s,x{s,z,0))yi{s,z)ds + o{£). 


So from (fT5|) . (13^ . and (|43|) we conclude that 
(44) 

x{t,z,£)= z+ Fo{s, x{s, z,0))ds 

Jo 

+£ / [D^Fo{s,x{s, z,0))yi{s, z) + Fi{s,x{s, z,0))]ds + o{£) 


= x{s, z,0) + £yi{t, z) + o{£). 

The last equality is a simple consequence of the computations made in Claim [TJ 
Indeed from fl26l) and Claim [T] if t^~^{z,0) <t< t^{z,0), then 

yi{t,z) = yi{t^~^{z,0),0)+ [ [D^Fo{s,x{s, z,0))yi{s, z)+Fi{s,x{s, z,0))]ds. 

Jt‘-^{z,0) 

From here, proceeding by induction on £, we obtain that 

yi{t,z)= / [D^Fo{s,x{s,z,0))yi{s,z) + Fi{s,x{s,z,0))]ds 
Jo 

This completes the proof of Claim [3] and, consequently, the proof of the lemma. 

□ 
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Lemma 3. Under the hypothesis of Theorem\^ there exists a compact subset Z 
of C with Z C Z° such that the solution x{t, z, 0) of the unperturbed differential 
system o is in the variable z for every z ^ Z. Moreover [dx/dz){t, z^U) = 
Y{t, z)Y (0, z)~^. The set Z is defined in the statement of LemmalE and Y is the 
fundamental matrix solution of (|8]) . 

Proof. It is easy to see that there exists a compact subset Z oi C such that 
Z <Z Z°. Given z ^ Z the solution of the unperturbed system ([7]) starting at 2 ; 
is given by ffTT]) by taking e = 0. Since C fl (9So = 0 and Z (Z C, there exists 
a neighborhood Z C oi z such that for every ( E the local flow of the 
unperturbed system ([7D starting at the point ( is given by a;^(t,C,0). We know 
that (t\z),x{f{z), z, 0)) G for i = 1, 2,..., — 1. Since is an open subset 

of S we conclude that there exist W Z neighborhoods of x{f{z), z, 0) in S for 
i = 1,2,..., Kz — 1. For i = Kz we have that x{U^{z), z, 0) = z, so we take = 
U^. Moreover, for each ( E W the locally flow of the unperturbed system ([7]) 
starting in ( is given by x^~^^{t, (, 0) for i = 1,2,..., Kz- Therefore we can choose 
a small neighborhood Uz Z such that for every E Uz, x(f(0,C,o) e 
for i = 1 , 2 , ..., Kz- Hence we conclude that for each z E C there exists a small 
neighborhood Uz Z C oi z such that the solution t i-A x(t, C, 0) can be written as 
(IT^ for every ( E Uz having the same number Kz of pieces. 

Let ipn(t,to,xo) be the solution of the differential equation x' = Ff(t,x) such 
that PnifoUoi^o) = From the results of the differential dependence of the 
solutions we conclude that each of these functions are of class in the variables 
it, to, xq). Indeed the function Ff is for i = 1, 2,..., Kz- From Claim [T] of the 
proof of Lemma [2] for i = 1,2,..., Kz the function U{C, e) is of class C^, for every 
C, eUz and e E [0, ^o]- 

From f[T5]) we have that 


(45) 


C,0) = 0,C) and 


for C E Uz and for i = 2, 3,..., Kz- So for i = 1 the function {t, () i-A x^(t, (, 0) = 

(Pni(i, 0, C) is C^. Moreover for 0 < t < t^((, 0) we have that -^—U, C, 0) = H (t, Q. 

oz 

Indeed from flTT)) we have that 


dx^ 

dz 



= D,Fo"Ht,x'(t,C,0)) —(t,C,0) 

dx^ 

= L>^Fo(t,a:(t,C,0)) —(t,C,0), 


( 46 ) 


Wt 
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for 0 < t < t^(C,0). So solving the linear differential equation (l46l) we have 
dx^ 

that the (t,C,0) is a fundamental matrix solution of system dHl) for 0 < t < 
az 

dx^ 

t^(C,0) and C G U^. Since ——(0,^,0) is the identity matrix, we conclude that 

yj Z 

f)x^ 

C, 0) = y(t, C)Y (0, z)-'^ for 0 < t < 0) and C e U^. 

We assume by induction hypothesis that the function C. •—t is 

dx^ 

for each t G and that for 0) <t < 0) the equality (, 0) = 

oz 

F(t,C)h^(0,z)-^ holds. 

From (1^ we have that, for i = i, x^{t, (, 0) = V5n^(t, 0), x^~^X~^((, 0), (, 

0)). So the the function ( i—)> x^(t, (, 0) is because from the induction hypothesis 
it is composition of functions. Now, we have 

m =/^.F;^(t,x^(t,c,o))^(t,C,o) 

dx^ 

= D^Fo{t, x{t, C, C, 0 ), 

for 0) < t < t^{C, 0). Solving this linear differential equation we get that 

^(t, C, 0 ) = r (t, OyX-\C 0 ), 0 ), c, o) 

= F(t,c)no,C)"\ 

for < t < and ( G Uz- The last equality comes from the 

induction hypothesis because 

— 0 ), C, 0 ) = 0 ), C, 0 ) = yX-\c, 0 ), C)>^( 0 , C)-'. 

The above induction proved that for every z & Z, x*(t, z,0) is a function 

dx^ 

in the second variable and ——(t, z,0) = Y(t, z)Y(0, z)~^, provided that < 

oz 

t < f. We conclude the proof of the lemma by observing that for z & Z and 
t G there exists I G {1,2, ...,k^} such that t^~^{z,0) <t< f{z,0), hence 
x{t, z,0) = x^{t, z,0). □ 

Lemma 4. Under the hypotheses of Theorem H] there exists a small parameter 
e G [0,£o] such that for every e G [0,F] the function z i-A- x{T,z,e) is locally 
Lipshchitz for z & Z. The parameter Eq is defined in the statement of Lemma\^ 
and the set Z is defined in the statement of Lemma\^ 
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Proof. From Lemma |3] we have that for each z & Z there exists a small neigh¬ 
borhood Uz C oi z such that the solution 1 1 —)■ x{t, (, 0) can be written as flT^ 
for every ( ^ Uz having the same number Kz of pieces. Therefore applying 
the result of the continuous dependence of the solutions on the parameters in 
each differentiable piece we conclude that for each z & Z there exists a small 
neighborhood lAz C Uz and a small parameter Sz G (0, £o] such that the solution 
t i-A x(tX,£) can be written as flT^ for every C, E Uz and for each e G ( 0 , 62 ] 
having the same number Kz of pieces. Since Z is a compact set we can choose 
e a minimal parameter of Sz for z E Z such that the above result holds for every 
e G [0,£]. 

Let 'il^nit,tQ,XQ,e) be the solution of the differential equation 
(47) x' = F^{t, x) = F^it, x) + sFfit, x) + x, e), 


such that 'ilJnito,to,xo,e) = xq. Clearly a^o^O) = ipnit,to,xo) which has 

been defined in Lemma El From the result of the continuous dependence of the 
solutions on the initial conditions we conclude that each of these functions are 
continuous in the variables (t, to, Xq). Indeed F" is a continuous function which is 
Lipschitz in the second variable for i = 1 , 2 , Kz- Moreover using the Gronwall 
Lemma (see, for instance, [28]) we conclude that 


(48) \\ilJn{t,si,si,e) - i>n{t,t2,Z2,e)\ \ < - t2 \ + e^'^\\xi - X2\\, 

for each t,si,S 2 E zi,Z 2 E Uz, and e G [0,£], where the constant L and M 
are defined in the proof of Lemma [2l From the flow properties of the solutions of 
system (H7|) we have that the equality 

(49) + S, to, X, e) = fjnit, to, ifnis + to, R, X, e),e) 


holds for every n = 1,2,..., N. 

Given ti,t 2 ,si,S 2 G §^, zi,Z 2 E Uz and e G [0, F] we can prove that the inequal¬ 
ity 


\Xn{ti,Si,Zi,e) - ^pn{t2,S2,Z2,e)\ \ < - t2\ + - Sa 

- Z2\\. 


(50) 
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holds for n = 1, 2,..., iV. Indeed, from (l48l) and (14^ 

Si,Zi,e) - 1pn{t2, S 2 ,Z 2 ,e)\\ = 

-1pn{tl,S2,1pn{t2 - k + S2, S2, Z2, s), e)\ \ < 

e^'^(^\\zi - knih -ti + S2,S2,Z2,e)\\ + M\si - S2I) = 

(^\\'kn(t 2 — ti + S2,t2 — ti + S2,Zi,e) — Ipnik ~ + S 2 ; - 52 , ^ 2 ; ^) | | 

+M|si - S2I) < 

(Iki ~ ^2!! + M\ti — ^2! + M\si — S2I). 


Again from 0131) we obtain 


(51) 


a;^(f,C,£) ='0ni(t,O,C,£) and 

x\t, C, 0) = kuiit, f ^(C, e), "^(C, e), C, e), e), 


for C, (z Uz and for i = 2,3,..., Kz- Thus from flSTl) for i = 1 the function 
x^{t,(^,e) = ipnj^(t,0,C). So from fl50|) we have that 


\\x^(ti,zi,e) - a:^(t2,^2,£)|| = HkmitiyO, zi,e) - kniihyO, Z 2 ,e)\\ 

< {\\zi - Z 2 \\ + M\ti - t 2 \), 

for every Zi,Z 2 G Uz, 0 < ti < k{zi,e), 0 < t 2 < k{z 2 ,s), and e G [0,e]. 

We assume by induction hypothesis that there exist constants Ai_x and 
such that 


\\x^ ^{ti,zi,e) - x^ ^{t2,Z2,e)\\< Ai_i\ti-t2\+Bi_i\\zi-Z2\\, 

for every ^ 1,^2 G Uz, k~‘^{zi,e) < k < k~^{zi,e), k~‘^{z 2 ,e) < ^2 < t^~^{z 2 ,e), 
and e G [0, e]. 

From fl5T]) we have, fori = ^, that x^{t,(,e) = 
e),e) for ( G Uz, k~^{(,e) < t < k{(,e) and e G [0,£]. So from induction 
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hypothesis we obtain that 
(52) 

\\x\ti,zi,e) - x^{t 2 ,Z 2 ,e)\\ = zi, e),e) 

-'lpnt{t2,t^~^{^2,£),X^~^{t^~^{z2,e),Z2, £),£)|| < 

Me^'^\ti - t2\ + Me^'^\t^~^{zi,e) - t^~^{z2,e)\ 

+e^'^\\x^~^{t^~^{zi,e),zi, e) - x^-^{t^-^{z 2 ,e), Z 2 , e:)|| < 

Me^^\ti - t2\ + e^^(M + Ae-i)\t‘^~^{zi,e) - t^~^{z2,e)\ + e^'^B^-iWzi - ^211 

for every Zi,Z 2 G t^~^{zi,e) < ti < t^{zi,e), t^~^{z 2 ,s) < t 2 < t^{z 2 ,s), and 
e e [0,£]. 

From Claim 1 of the proof of Lemma [2] we have that t^~^{z, e) is a function, 
then there exists a constant 5 > 0 such that \t^~^{zi,e) —t^~^[z 2 ,e)\ < 5||zi — ^211 
for every e G [0,F]. Going back to the inequality fl52ll we get 

\\x\ti,zi,e) - x\t2,Z2,e)\\ < Ai\ti -t2\ + Be\\zi -2r2||, 

for every 2 : 1 , ^2 G Uz, t^~"^{zi,e) < A < t^{zi,e), t^~^{z 2 ,s) < t 2 < t^{z 2 ,s), and 
e G [0,F], where Ai = and Bi = {6{M + Ai^i) + 

We conclude the proof of the lemma by observing that a;(T, z, e) = x^^(T, z, e) 
which, from the above induction, is locally Lipschitz in the variable z. □ 

Lemma 5. Under the hypothesis of Theorem\M the solution x{t,z,e) of the un¬ 
perturbed differential system o is in the variable z for every z & Z. More¬ 
over {dx/dz){t,z,Gi) = Y{t, z)Y{0, z)~^. The set Z is defined in the statement of 
Lemma 0 and Y is the fundamental matrix solution of ()5|) . 

Proof. Assuming the hypothesis (hi) instead {HI) we can prove analogously to 
Claim [T] in the proof of Lemma [2] that given z & Z the function P{z,e) for 
i = 0,1, 2, • • • ,Kz is of class for every C neighborhood Uz <Z C oi z and 
£ G [0, £ 0 ]- Then the proof of the lemma follows analogous the proof of Lemma [3] 
but considering the functions 'fin{t,to,xo,e) defined in Lemma 01 □ 

The next two lemmas are versions of the so called Lyapunov-Schmidt reduction 
for finite dimensional function (see for instance m) and its proof can be found 
in [9] and [8] , respectively. The hrst lemma will be used for proving Theorem 
and the second one will be used for proving Theorem [Bl 

Lemma 6. Let P : ^ be a function, and let Q : Efi x [0,eo] ^ be 

a continuous functions which is locally Lipschitz in the first variable, and define 
f : X [0,£o] f{z,e) = P{z) + eQ{z,e). We assume that there exists 

an open and bounded subset 1/ C R^ with k < n and a function : V ^ R'^”^ 
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such that P vanishes on the set Z = {zq, = (a,/3o(Q!)) : a G V} and that for 
any a the matrix DP{za) has in its upper right corner the null k x {d — k) 
matrix and in the lower corner the {d — k)x{d — k) matrix A„ with det(Aa) ^ 0. 
For any a E V we define /i(a) = 7rQ{za,0). Thus if fi{a) 7 ^ 0 for all a G dV 
and d{fi, V, 0) 0, then there exists £1 > 0 sufficiently small such that for each 

e G (0,£i] there exists at least one z^ G with F{zs,e) = 0 and dis{ze,Z) — )■ 0 
as £ —)■ 0 . 

Lemma 7. Let P ^ and Q ■.'Mfix [0, £ 0 ] functions, and define 

f : 'Mfi X [0,£o] —^ as f{z,e) = P{z) + sQ{z,e). We assume that there exists 

an open and hounded subset C with k < n and a function fio : V 
such that P vanishes on the set Z = {za = («,/9o(a)) : a G V} and that for any 
a eV the matrix DP{za) has in its upper right corner the null kx{d — k) matrix 
and in the lower corner the {d — k) x {d — k) matrix A„ with det(Act) 7 ^ 0. For 
any a E V we define fi{a) = 7iQ{za, 0). Thus if there exists a E V with fi{a) 7 ^ 0 
and det(/'(a)) 7 ^ 0 , then there exists as such that f{za^,e) = 0 and Za^ Za as 
e ^ d. 

Now we are ready to prove our main results. 

Proof of Theorem\M We consider the function f : Z x [0,£o] given by 

(53) f{z,e) = x{T,z,e) - z. 

Its differentiability comes from Lemma [HI Clearly system ([ 6 ]) for e = e E [0,£o] 
has a periodic solution passing through 'z E C li and only if f{z,e) = 0. 

From Lemma [2] we have that 

(54) x{t, z, e) = x{t, z, 0) + eyfit, z) + o{e). 

Taking P{z) = x(t, z,0) — z and Q{z, e) = yfit, z) + o(e), thus f{z, e) = P{z) + 
eQ{z,e). Moreover from Lemma [3] P(. 2 ) is a function, and from Lemma |4| 
Q{z,e) is a continuous function which is locally Lipschitz in the hrst variable 
because Q{z, e) = {x{T, z, e) — x(T, z, 0))/e. 

In order to apply Lemma [ 6 ] to function fl53D we compute 

P{Za) = X{T, Za, 0)- Za= 0, 

and 

3 P St 

^(^«)= Y^{T,Za,0)-Id 

= Ya{T)Ya{0)-^ - Id. 

So from hypothesis (H) the function P vanishes on the set Z and from hypothesis 
{H2) for any a E V the matrix DP{za) has in its upper right corner the null 
k X {d — k) matrix and in the lower corner the {d — k) x {d — k) matrix Aq, 
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with det(A„) 7 ^ 0. Since 7 rQ(a,/So(a)) = -KyiiT^Za) = /i(a), so the proof follows 
applying Lemma Ei □ 

Proof of Theorem O The proof is analogous to the proof of Theorem |A] applying 
Lemma |5] instead of Lemmas |3] and IU and applying Lemma [7] instead of Lemma 

E □ 


3. Proof of Proposition □ 

Proof of Proposition E Proceeding with the change of variables (n, n, w) = (r cos 9, 
r sin 6 ', z) and taking 9 as the new time by doing r' = r/9 and z'zl9 we obtain 


(56) 


{0, z) + e G^{9,r, z) + 0{e‘^) if 0 < 0 < tt, 
{0, z) + e G~ {9,r, z) + if 7i<9<27i, 


where G^ = [G^,Gf), and 


Gf = b^r cos^ 9 + (a^ + dfz + {bf + c^)r sin 9) cos 9 
+ [af + dfz + cfr sin 9) sin 9, 

Gf = - [r{af + dfz) — bfrz cos^ 9 + (Gfr"^ + (a]*^ + dfz)z + cfr sin 9) sin 9 
{bfr'^ — {af + df)z + {bf — cf)rzsm9) cos 9) . 

Here the prime denotes the derivative with respect to 9. 

For system (1551) we have that D = {{r, z) : r > 0, z G M} and T = 27r. We 
note that S = {(0,r) : r > 0} U {(vr, r) : r > 0} U {(27r,r) : r > 0}, thus taking 
h{9,r,z) = 9{9 —'k){9 — 27r) it follows that S = h“^(0). 

In what follows we shall study the elements of hypothesis {H) of Theorem lAl 
For e = 0 the solution x(0,r, z,0) of system (l55|) such that x{0,r,z,0) = {r,z) 
is given by x{9,r,z,0) = {r,e^z). Taking P = {rGM: ri < a < r 2 } with 
ri > 0 arbitrarily small and r 2 > ri arbitrarily large, and /3o = 0 we have that 
the solution Xa{9) = (a,0) is constant for every a G P, particularly 27r-periodic. 
In this case the manifold Z of periodic solution of the system fl55|) when e = 0 
is given hj Z = {(a,0) : ri < a < r 2 }, and Sq = D. Since Z C Sq it follows 
that Z n 9So = 0. Moreover computing the crossing region of system fl55|) for 
£ > 0 sufficiently small we conclude that = S, so we obtain that Zq fl S C 
Therefore hypothesis (FT) hods for system (155|) . 
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Hypothesis {HI) of Theorem lAl clearly holds for system (j55l) . To verify hy¬ 
pothesis {H2) we take 


3t 

Y{0,r,z) = —= 


1 0 
0 e® 


as the fundamental matrix solution of system ([8]) in the case of system (l55|) . So 


y„(27r)y„(0)-^ -Id = Y(27r, a, 0)^(0, a, 0)"^ - Id = 


0 0 
0 e^^-l 


Since Aq = — 1 7 ^ 0 for every a G H it follows that hypothesis (H2) holds for 

system fl^ . 

Now if (0, r, z) G S, then 6 G {0, vr}. On the other hand Vh(0, r, z) = (27r^, 0, 0) 
and Vh{'K,r,z) = (— 7 r^, 0 , 0 ) for every {r,z) G D. So (Vh( 6 *, r, 2 ;), (0, n)) = 0 for 
every 9 G { 0 , 7 r}, {r,z) G D, and v G which means that for any n G we 
have that (0, v) G T (0 for every 6 G {0, vr} and (r, z) G D. In short hypothesis 

{H3) holds for system ([55]). 

Using an algebraic manipulator as Mathematica or Maple we compute 
/i(«) = ^{bt + K + ct + C2)a + 2 [at - a^) . 

From hypotheses + h'^ + C 2 + c^) (a^ — >0, thus 

_ 4 (®2 ~ ^2 ) 

vr ( 6 ]^ + h'^ + 4 + C 2 ) 

is a solutions of the equation fi{a) = 0 such that f[{a) 7 ^ 0. From Remark [2] it is 
a sufficient condition to guarantee the existence of a small neighborhood kF C U 
of a such that d{fi,W,0) 7 ^ 0. Since /i is linear, it is clear that fi{a) 7 ^ 0 for 
every a G dW. Therefore hypothesis (HA) of Theorem lAl holds for system fl55l) . 

Now the proof of the proposition follows directly by applying Theorems IXI and 

El □ 
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